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ABSTRACT 

Aims. We want to study the temporal and spectral behaviour of HU Aqr in the X-ray domain during different accretion 
states. 

Methods. We obtained spectra and light curves from four different XMM-Newton pointings covering intermediate and 
low states. The X-ray observations were accompanied with high time resolution photometry obtained with the Optima 
and ULTRACAM instruments. 

Results. On two occasions in May 2002 and 2003 HU Aqr was found in an intermediate state with the accretion rate 
reduced by a factor of 50 compared to earlier high state measurements. X-ray spectra in the intermediate state can be 
described by a model containing a blackbody component and hot thermal plasma. Contrary to the high state the ratio 
between soft and hard X-ray flux is nearly balanced. In agreement with previous measurements we observed a migration 
of the accretion spot and stream towards the line connecting both stars. The brightness of HU Aqr was further reduced 
by a factor of 80 during two low states in October 2003 and May 2005, where it was detected at a luminosity of only 
Lx = 4.7 x 10 28 erg s~ : . This luminosity would fit well with an active coronal emitter, but the relatively high plasma 
temperatures of 3.5 and 2.0 keV are more compatible with residual accretion. We updated the eclipse ephemeris of 
HU Aqr based on the eclipse egress of the accretion spot measured in various wavelength bands. The (O — C)-diagram 
of the observed accretion spot eclipse timings reveals complex deviations from a linear trend, which can be explained 
by a constant or cyclic period change or a combination thereof. The quadratic term implies a period decrease at a rate 
ofPorb = -7..-11 x 10" 12 ss" 1 . In case the observed period change reflects a true angular momentum loss, this would 
be a factor of 30 larger than given by gravitational radiation. 

Key words. X-rays: stars - stars: cataclysmic variables - accretion - stars: magnetic fields - AM Her systems - stars: 
individual: HU Aqr 



1. Introduction 

Polars or AM Herculis stars are important objects to un- 
derstand accretion in the presence of a strong mag netic 
field. In these Cataclysmic Variables (Warner 199§) the 
field of the white dwarf primary is sufficient to lock its 
rotation to the binary orbit, and to prevent the formation 
of an accretion disc. To date, we know 15 eclipsing sys- 
tems of this kind, which additionally provide constraints 
on the viewing geometry and allow the different emission 
sites to be disentangled. Among those, HU Aqr is one the of 
the brightest and has become one of the most comprehen- 
sively s tudied systems. Since its discovery in the ROSAT 
survey (jSchwope et al.l 119931 ) it has been studied in vari- 
ous wavelength bands, and was the first system in which 
the ballistic stream w as detected via Doppler tomography 
(jSchwope et al.|[l997l ). It features a rather simple accretion 
geometry with one dominating accretion spot that under- 
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goes self-eclipses. It was the target of a long-term moni- 
toring programme w ith the ROSAT and EUVE satellites 
(jSchwope et al.ll200l[ ). Part of these observations were ob- 
tained in a high state that revealed a strong soft X-ray 
excess and could be used to resolve the eclipse egress, thus 
directly constraining the diameter of the soft X-ray emit- 
ting area to ~ 450km. Modelling of the X-ray eclipse dur- 
ing the high state located the accretion region at azimuth 
ijj = 46° and colatitude (3 = 31° and required a substantial 
amount of vertical extent. During the monitoring HU Aqr 
displayed large variations of the accretion rate, accompa- 
nied by a correlated migration of the accretion spot and 
the coupling region. For the reduced (intermediate) state 
of accretion both accretion spot and stream were shifted 
towards the line connecting both stars in agreement with 
an increased threading radius given by the magnetic/ram 
pressure balance. All of the intermediate state observations 
where obtained with the ROSAT-HRI and EUVE only, 
thus lacking spectral information. In this paper we present 
a series of XMM pointings that continue the monitoring 
through intermediate and low brightness states of the sys- 
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Table 1. Log of the XMM X-ray observations of HU Aqr 



Date 


Detector^ ' 


J Exp 

(ksec J 


(2) 

Mean rate^ ' 
(cts/s) 


2002 May 16-17 


PN-sw 


36.4 


1.11 




MOSl-sw 


37.2 


0.19 




MOS2-sw 


37.2 


0.26 


2003 May 19 


PN-ti 


18.5 


0.94 




MOSl-ti 


20.2 


0.21 




MOS2-ti 


20.0 


0.18 


2003 Oct 24-25 


PN-sw 


11.4 


0.011 




MOSl-sw 


11.6 






MOS2-sw 


11.6 




2005 May 16 


PN-ti 


18.6 






MOSl-sw 


19.4 


0.0025 




MOS2-sw 


19.4 


0.0022 



' 1 ' The abbreviations after the PN/MOS detector indicate the 
instrument mode: sw = small window, ti = timing 
^ In the bright phase interval <f> = 0.62 — 0.91 and 4> = 
0.05 - 0.18. Energy range used 0.18-10 keV and 0.3-10 keV for 
imaging and timing modes, respectively. Count rates have been 
corrected for the 65% encircled energy fraction corresponding to 
the aperture used. 

tern. Preliminary res ults of this stud y were presented in 
ISchwope et ail (|2004f i and IVogel et"aH (|2008l ). 

2. Observations and reduction 

Between 2002 and 2005 HU Aqr was the target of four ded- 
icated XMM pointings. Different instrument configurations 
of the EPIC cameras (Tabled]) were used to best match the 
brightness level of the source. 

For the MOS and PN data taken in image mode source 
photons have been extracted within a circular aperture of 
15". The resulting count rates given in the paper have been 
corrected for the 65% encircled energy fraction correspond- 
ing to this radius. The background level has been deter- 
mined from a larger source-free region close to the source. 
For the timing mode observation the level of background 
radiation was comparable to the source flux resulting in 
much poorer data quality. The PN-data additionally suf- 
fered from soft proton flares, so here all photons with ener- 
gies lower than 0.3 keV had to be omitted. Apart from this 
particular data set, spectral analysis was applied to the 
energy range between 0.18-10 keV including good events 
(flag = 0) identified with patterns 0-12 and 0-4 for the 
MOS and PN-CCD, respectively. We used version 7.0.0 of 
the XMMSAS package. 

During each pointing the Optical Monitor (OM) was 
operated with the filter UVM2 using the fast mode. These 
data were reduced using the task OMFCHAIN. For conversion 
of the count rates to flux units we applied the factors given 
in the XMM calibration file OM.COLORTRANS. 

We also obtained contemporary optical data with 
the fast-timing p hotometric devices ULTRACAM 
(jDhillon et all l2007f i and Optima (|Kanbach et al.ll2003h . 
These instruments were operated at the 4.2-m William 
Herschel Telescope (WHT) on La Palma, the VLT 
on Paranal, and the 1.3 m telescope of the Skinakas 
Observatory, Crete. A detailed description of these optical 
and the OM UV observations will be given in a different 
paper (Vogel et al., in prep.). 
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Fig. 1. Upper panel: Residuals of the timings of accre- 
tion spot egre s s with respect to the linear ephemeris of 
ISchwope et alj (|2001l ) together with the best-fit quadratic 
ephemeris derived in this paper (solid line). Middle panel: 
Residuals with respect to our new linear ephemeris to- 
gether with the best-fit linear+sinusoidal ephemeris (solid 
line). Lower panel: Residuals with respect to the linear and 
quadratic terms of the combined quadratic and sinusoidal 
ephemeris. Overplotted as a solid line is the sinusoidal term 
of this ephemeris. 

3. Results 

3.1. An updated accretion spot ephemeris 

We have updated the eclipse ephemeris of HU Aqr using 
newly available timings of the spot egress measured in the 
optical, UV and X-rays (Ta ble d and the older ROSAT 
and EUVE measurements of ISchwope et al.l (|200lh . As the 
timestamp of egress/ingress in the optical we used the mo- 
ments of half intensity. For the photon counting data (UV 
and X-ray) this datum was computed from the mean of 
the arrival times of the first three photons after the eclipse. 
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Table 2. Eclipse egress ephemerides of HU Aqr 



Linear ephemeris: 

BJED = T + P ■ E 

T = 2 449 9102.9200839 (±6) 
xl = 345.9, v = 101 

Quadratic ephemeris: 

BJED = To + P ■ E + Q ■ E 2 

To = 2 449 102.9200164 (±8) 
Q = (-3.17 ±0.02) x 10~ 13 
Xl = 40.5, v = 100 



Po = 0.08682041087 (±3) d 



Po = 0.08682042416 (±8) d 



Sinusoidal ephemeris: 

BJED = T + Po • E + A- sin [2n(E - 

T = 2 449 102.9201788 (±15) d P 
A = (125 ± 1) x 10~ 6 d B = 

C = (554 ± 5) x 10 2 cycles 
xl = 30.4, v = 98 



B)/C] 

= 0.08682040612 (±6) d 
(89 ± 4) x 10 3 cycles 



Quadratic + sinusoidal ephemeris: 

BJED = T ± Po ■ E ± Q ■ E 2 ± A- sin [2ty(E - B)/C] 

To = 2 449 102.9198961 (±25) Po = 0.0868204353 (±2) d 



Q = (-4.77 ±0.04) x 10~ J 
B = (-4190 ± 150) cycles 
X2 = H, ^ = 97 



A 
C 



(101 ± 2) x 10" 6 d 
(291 ± 1.5) x 10 2 cycles 



Excluded were the low state OM light curves, which are of 
low counting statistic and likely dominated by the white 
dwarf photosphere. All timings have been converted to the 
solar system's barycenter and are in the time frame of 
ephemeris time which accounts for leap seconds. 

In Fig. [1] we compare the observed egre ss timings with 
the pr edictions of the linear ephemeris of ISchwope et alj 
(200lJ), which reveals conspicuous deviations of up to 45 
sec for the latest observations. We have updated the linear 
ephemeris using a least-squares fit weighted by the inverse 
of the squares of the eclipse timings. The best-fit parame- 
ters together with their formal l-a errors are given in Table 
[2j The fit is very poor, showing large positive and negative 
residuals. 

A better match to the data can be achieved by introduc- 
ing either a quadratic or a sinusoidal term. The first would 
indicate a secular variation of the orbital period, while the 
latter indicates a cyclical period change. The quadratic 
term is large and would imply an orbital period decrease 
of P or b = —7.3 x 10~ 12 s s~l, which is larger than what is 
expected from angular momentum loss due to gravitational 
radiation. The cycle length of the sinusoidal modulation 
is 13 yr, close to the time-span of our observations. The 
quadratic and sinusoidal ephemerides are shown as solid 
lines in the upper and middel panels of Fig.[TJ Both models 
feature a large Xv of 40.5 and 30.4 emphasizing that the 
period variation is more complex than these simple fits. Of 
particular note is a rapid increase in the (O — C)-values of 
10 sec/yr between 1998 and 2000, and a subsequent strong 
period decrease in the following years lasting until 2004. 
This could be indicative of a cyclical period modulation 
on a shorter time-scale of just a few years. As a possible 
model we assumed a combination of a quadratic plus sinu- 
soidal term, which account for a long-term period change 
and cyclical modulation simultaneously. The result of this 
fit is displayed in the lower panel of Fig. [TJ where the resid- 
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Fig. 2. upper panel: Phase-averaged EPIC PN X-ray light 
curve of HU Aqr observed with XMM-Newton in May 2002. 
The black and gray line mark the light curve in the hard 
(0.5-10 keV) and soft (0.18-0.35 keV) energy bands, lower 
panel: Variation of the hardness ratio ((H-S)/(H+S)) as a 
function of orbital phase. 



uals of the linear and quadratic part of the ephemeris are 
shown together with the fitted sinusoidal term. The com- 
bined fit yields a P or b of —1.1 x 10~ n s s _ l, slightly larger 
compared to the simple quadratic ephemeris. The cyclical 
variation has a period of 6.9 yr and an semi-amplitude of 
~ 9 sec. Our data so far faithfully revealed two consecutive 
cycles of this putative modulation lending some support for 
its existence. On the other hand, the very latest data points 
are already deviating by 10 sec suggesting that the derived 
period decrease might be to large. 

With a Xv of 11.5 the fit is clearly better than the single 
quadratic or sinuso idal models. U sing an F-test following 
the prescription of iPringld (|1975[ ) the significance of the 
quadratic plus sinusoidal ephemeris compared to the linear 
plus sinusoidal ephemeris is 93.9 per cent with P(l,97) = 
170. While this model is statistically favoured, it is too early 
to say what the future behavior of the period wanderings in 
HU Aqr will be. However, the eclipse timings to be taken 
in the next few years should clearly discriminate between 
the different possible solutions and confirm the existence of 
a large quadratic term. 

3.2. X-ray and UV light curve morphology 

In Fig. [5] and Fig. [3] we show the EPIC-PN light curves dur- 
ing the intermediate states observed in May 2002 and May 
2003. For the data from May 2002 we extracted photons in 
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Fig. 3. Phase-averaged EPIC PN X-ray light curve of 
HU Aqr in the 0.3-10 keV range observed with XMM- 
Newton in May 2003. 
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Fig. 4. Phase-averaged OM light curves observed in the 
UVM2 filter at 2200 A. The two upper curves are from 
the intermediate states of May 2002 (black) and May 2003 
(grey), while the two lower graphs correspond to the low 
states of October 2003 (dots) and May 2005 (histogram). 
The flux unit is 10~ 15 erg cm -2 s -1 A -1 . 



two spectral bands ranging from 0.18-0.35 and 0.5-10 keV, 
which match the two spectral components (see below). A 
corresponding hardness ratio has been calculated with the 
two bands related as (hard-soft)/ (hard+soft). 

As in previous X-ray observations the light curves are 
dominated by a single bright phase interval of a self- 
eclipsing pole. The centre of the bright phase of the main 
accretion region was at <f> = 0.912 and cf> = 0.915 for the 
May 2002 and 2003 observation, respectively. This corre- 
sponds to an azimuthal position of tp ~ 31°. There is no 
obvious difference between the location of the soft and hard 
X-ray emitting region. Compared to previous observations 



these new measurements show the main accretion region 
most closely shifted towards the line connecting both stars. 
Such a relation between position and instantaneous accre- 
tion rate is expected for an intermediate state from the 
gener al trend seen in some older ROSAT/EUVE observa- 
tions (jSchwope et al.ll200lh . 

The light curve from May 2002 reveals the presence of 
narrow and broad absorption dips in the soft X-ray light 
curve. For the observation of May 2003 the counting statis- 
tic in the soft band is not sufficient for a clear detection of 
either dip. In May 2002 the narrow absorption feature was 
centred at phase <j> — 0.941 with the dip ingress starting 
at (j> = 0.908. These values are among the closest to the 
eclipse observed, indicating a migration of the stagnation 
point towards the line connecting both stars. A broad ab- 
sorption trough is centred around <f> = 0.73 very similar to 
what was previously observed in several ROSAT and EUVE 
pointings. 

There is a substantial hardness ratio variation during 
both the narrow and the broad absorption dip. The nar- 
row dip is not total even in the softest X-ray band. A fit 
to the dip spectrum with the temperature of the black- 
body components fixed yields a cold absorption column of 
(2.7± 0.2) x 10 20 cm" 2 . This is a factor often lower then the 
(2 - 3) x 10 21 cixr 2 measured during the ROSAT high-state 
observation from October 1993. 

Contrary to the bright-state ROSAT observation from 
October 1993, eclipse ingress and egress profiles are un- 
resolved for the XMM observations at a timescale shorter 
than ~ 2 sec due to the low count rate of the data. 

Fig. 2] shows the phase folded light curves of HU Aqr 
taken with the Optical Monitor and the UVM2 filter. 
The source was in two different brightness states, be- 
ing brighter by a factor of two during the intermediate 
state observed in May 2002 and May 2003. During both 
low states HU Aqr was at the same brightness level, be- 
tween 0.8 x 10 _15 erg cm -2 s -1 A -1 in the faint phase and 
2 x 10 _15 erg cm~ 2 s _1 A -1 at orbital maximum. At all 
four epochs we observe a smooth sine-like modulation rem- 
i niscent of a heated c ap on the surface of the white dwarf 
(Gansick e et al.l [20061) . The level of modulation is rather 
similar for all brightness states. This suggests that the 
amount of heating does not strongly depend on the instan- 
taneous irradiation. A detailed modelling of this feature will 
be presented in Vogel et al., in prep. 

In the intermediate brightness state the pre-eclipse UV 
flux is depressed by 30% starting at phase 0.92 (see inset of 
Fig. 2]) . A plausible explanation for this dip is absorption 
by the accretion stream also seen in soft X-rays. The addi- 
tional brightness compared to the low state is seen in the 
bright and faint phase, i.e. it does not come from the self- 
eclipsing heated cap. A possible source of this emission is 
the accretion stream, which is also indicated by a prolonged 
ingress after eclipse of the white dwarf. 

3.3. X-ray spectroscopy 

X-ray spectra in the intermediate state have been extracted 
using the bright phase intervals from <fr = 0.62 — 0.91 and 
(f> = 0.05 — 0.18, thus excluding the absorption dip and the 
eclipse. The spectra of the PN and MOS detectors have 
been fitted together and can be described with a two com- 
ponent model including a soft blackbody and hot thermal 
plasma. For the latter both a mekal or a bremstrahlung 
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Fig. 5. EPIC-PN bright phase spectrum of HU Aqr taken 
in May 2002 together with a fit containing a blackbody and 
mekal model modified by ionised absorption and cold re- 
flection. Center: Residuals with respect to the fit. Bottom: 
Relative contribution of the ionised absorber (dashed line) 
and reflection component (solid line) to the total X-ray 
spectrum. 



model fits relatively well with a x 2 of 1.46 and 1.49. The 
observation from May 2002 is sufficient to resolve individ- 
ual line features. Above 6 keV all three lines of fluorescent, 
H-like and He-like iron are present and do contribute to 
the x 2 - We also note weak residuals seen at 0.5 keV from 
a blend of line emission, probably dominated by the O vn 
and Ovili lines at 21.8 and 19 A. The ratio between the 
H-like and He-like iron lines indicates plasma with tem- 
peratures between 5-15 keV, while the continuum slope at 
energies larger 5 keV forces the fit to unrealistically high 
temperatures of fcTb rem s = 115 keV and fcT me k = 98 keV. 
Introducing a continuous temperature CMEKAL model does 
not help to solve the discrepancy between the temperature 
indicators. 

A way to reconcile the spectral slope with more realistic 
plasma temperature is by adopting an ionised absorber as 
imple mented in XSPEC by the ABSORI model ()Done et al.l 
|1992| ) . Introduction of the ABSORI model significantly im- 
proves the fit to xl — 1-38. Using an F-test this model 
is better than the model without warm absorption at the 
75% significance level. Comparison with fits to the PN- 
data alone yields a \1 of 1.12, showing that the remaining 
discrepancy is mainly due to the uncertainty in the cross- 
calibration of the different detectors. As parameters of the 



ABSORI model we fixed the power law index of the incident 
radiation to 1.25 and the temperature of the ionised mate- 
rial to 30000 K. The remaining best fit parameters then are 
an absorption column of 6.5 x 10 21 cm -2 and an ionisation 
parameter £ = 110 erg cm sec -1 . For the observation of 
May 2003 these parameters were fixed at the values from 
2002 given the low quality of the data. The complex ab- 
sorption drastically decreases the best-fit temperatures of 
the mekal model to 19 and 15 keV for the observation in 
2002 and 2003, respectively. 

The incidence of the fluorescent iron line at 6.4 
keV reveals that there is also reflection from neutral 
material at the surface of the white dwarf. This has 
been taken into account by adding the reflect model 
(|Maed ziarz fc ZdziarskJ [l995) to our fitting. The new 
model component only marginally improves the x 2 , but 
lowers the temperatures of the mekal component to 15 
and 13 keV for the observation in 2002 and 2003, respec- 
tively. The only free paramter is the scaling factor of the 
reflected component, which for the fit of the 2002 data is 
R = 1.66. This value was fixed for the fit to the much poorer 
data from 2003. For completeness we also provide a fit using 
only the reflection model without an ionised absorber. For 
the 2002 observation this fit has a significantly higher % 2 
compared to the fit including warm absorption, while the 
difference between both models is negligible for the data 
of May 2003. Using this model the best-fit temperatures 
of the MEKAL component change to 35 and 22 keV for the 
observation in 2002 and 2003, respectively. 

The blackbody temperatures for the observations of 
May 2002 and May 2003 differ and also depend on the 
details of the absorption and reflection models used. The 
best-fit values for our finally adopted model including warm 
absorber and reflection are fcTbb = 34 eV and fcTbb = 25 
eV, differing only slightly above the 2a-level. Those values 
fall right into the range of temperatur es measured for po- 
lars with a dominant soft component (iRamsav fc Cropper! 
120031 : IRamsav et al.l [20041 ) . The decrease of the tempera- 
ture in May 2003 is accompanied by an increase of the 
unabsorbed, bolometric soft X-ray flux by a factor of 6. 
At the same time the bolometric hard X-ray flux decreased 
by 10% . A much more significant change (at the 3cr level) 
of the blackbody temperature is observed if one compares 
the 34 eV of the intermediate state in 2002 with the 22 eV 
measured during the previous high state in October 1993 
(jSchwope et al.ll2001h . 

We also note the change of the hydrogen column density 
compared to the high state ROSAT observation. Both inter- 
mediate state observations require at maximum a density 
of 6 x 10 19 cm -2 , whereas the ROSAT pointing is consistent 
with a larger column of 2 x 10 20 cm~ 2 . The large variabil- 
ity of the absorption column point toward an intra-binary 
origin. 

Using the values from the fit including absorption and 
reflection, and taking it and 2ir for the scaling of the 
blackbody and thermal emission, the unabsorbed flux ra- 
tio /boi, soft //boi, hard were 1.3 and 8.4 for the May 2002 and 
May 2003 intermediate states. Taking into account the un- 
certainty due to the unknown emission characteristic and 
the neglect of the cyclotron component, the value from May 
2002 is relatively close to th e ratio of ~ 0.5 of the stan- 
dard accretion s hock m odel of lLamb fc Masters! |l979) and 
King fc Lasotal |l979). Both values mark a dramatic de- 
crease of the soft X-ray excess compared to the high state of 
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Fig. 6. Light curves of the simultaneous XMM- 
MOS/ULTRACAM observation of HU Aqr on May 16, 
2005. The relative flux units of the optic al light curves 
are w ith respect to the comparison star of ISchwope et all 
(1993). The blue, green and red lines denote the light curves 
in the u', g' and r' filters. 



October 1993, where the soft-to-hard ratio of the bolomet- 
ric fluxes was ~ 800. Note that this change is almost solely 
due to the decrease of the soft X-ray component, whereas 
the hard X-ray flux remained nearly constant from the high 
to the intermediate state. 

The total luminosity in X-rays Lx = (ftfbb + 
27r/therm)-D 2 , assuming a distance of 180 pc, was 2.5 x 10 31 
and 9.3 x 10 31 erg s" 1 in May 2002 and May 2003. The 
inferred accretion rate is between (0.8 — 3.1) x 10 _11 M Q /yr 
, a factor of 20 lower compared to the high state. 

3.4. Deep low-state observations 

In October 2003 HU Aqr was found in a deep X-ray low- 
state. An optical light curve taken on October 21, 2003 
three days in advance of the X-ray observations showed 
HU Aqr with a bright phase magnitude of only R = 17? 1 4 
. The source was detected with the PN detector at a rate 
of only 0.011 ± 0.003 cts/sec. The spectrum of the ~ 90 
source photons can be fitted with a thermal plasma mekal 
model with a xt °f 0-5- A fit using Cash statistics led 
to similar results. No dedicated soft blackbody component 
is required for the low state of HU Aqr. The fit temper- 
ature of 3.5 ± 3 keV is not very well constrained. The 
unabsorbed bolometric flux of the observed spectrum is 
4.8 x 10~ 14 erg cm' 2 s _1 yielding an X-ray luminosity of 
4.6 x 10 28 erg s _1 for a distance of 180 pc. There is weak in- 
dication for a variability pattern, where the photons cluster 
in the X-ray bright phase. 



A second deep low state was observed in May 2005. This 
observations was accompanied by a series of ULTRACAM 
light curves taken on consecutive nights between May 14 
and 17, 2005. In the g' and r'-band residual accretion ac- 
tivity was evident by cyclotron beaming and flaring. The 
w'-band, on the otherhand, was clearly dominated by sine- 
like variations reminiscent of a heated white dwarf cap. 
Starting from a bright phase magnitude of r' = 1T?5 seen 
on May 14, 2005 the level of additional accretion activity 
increased from night to night. For the simultaneous ob- 
servations on May 16, 2005 the maximum brightness was 
at r' = 17^3. HU Aqr was detected in both MOS detec- 
tors at mean rates of 0.0022 and 0.0025 cts/sec. The X-ray 
light curve (Fig. [5]) reveals clear variability, which was only 
partially correlated with the visibility of the main accre- 
tion region. For instance the source appears X-ray bright 
during one faint phase interval. As for the first low state 
observation the spectrum can be sufficiently fitted by a sin- 
gle temperature thermal plasma model. The fit parameters 
imply a negligible amount of absorption and a temperature 
of fcT mo k = 2 ± 1 keV. The source was at a similar bright- 
ness level compared to the October 2003 low state with a 
bolometric flux of 4.3 x 10 _14 erg cm~ 2 s _1 . 

The observed X-ray emission c an be compared wit h 
that of single, active late type stars (jPizzolato et alj |2003). 
Assuming a bolometric l uminosity of 3 . 9 x 1 31 erg s _1 
for the M4.5 secondary (jLeggett et al.l Il996f ) results in 
Lx/Lboi = 1.2 x 10~ 3 . This value is compatible with the 
saturated X-ray regime observed for the fastest rotators, 
although it falls into the top end being a factor of two 
higher than the mean. An argument against a coronal ori- 
gin are the relatively high plasma temperatures, that lie 
above the peak of log(T(K) ) = 7 in the temperature d istri- 
bution of coronal emitters ()Sanz-Forcada et all 120031 ) and 
are also higher than the tempera tures of ~ 0.6keV ob- 
served in low accretion rate polars (jVogel et al.ll2007f ). The 
plasma temperatures measured in the low-state of HU Aqr 
are better matched with an origin in a low density accre- 
tion c olumn that predominantly cools via cyclotron radi- 
ation (|Woelk fc Beuerm ann 1992). Signs of residual accre- 
tion are obviously present in the r' and g' light curves in 
May 2005, where cyclotron beaming and flaring is observed. 
The problem with this interpretation is that one would ex- 
pect a clear on/off variability pattern in X-rays due to the 
self-eclipse of the accretion region, which is not observed at 
least for the low state in May 2005. 

4. Discussion 

One of the major results of our study is the dependence 
of the soft X-ray excess on the X-ray brightness level and 
thus the instantaneous accretion rate. Unlike the high state 
when a strong soft X-ray excess was observed, the system 
reverted to a more balanced flux ratio in its intermediate 
accretion state. A similar result has been reported for a 
handful of systems in the X MM snap-shot survey of polars 
()Ramsav fc Cropper! I2004T ). It can be understood in terms 
of a density change in the coupling region, and is therefore 
an indirect confirm ation of the blobby accretion scenario 
([Frank et alJ ll.988). Our observation indicates a tempera- 
ture increase of the blackbody component from 22 eV to 
34 eV as the system changed from the high to the inter- 
mediate state. There are only few observations of polars 
which constrain the physical parameters that determine the 
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Table 3. Fit results for the combined PN and MOS X-ray spectra of HU Aqr 
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Model JVh fcTbb fcT thcrm Fsoft.bol -Fhard.bol Xred 

(10 19 cm" 2 ) (eV) (keV) (erg cm^s" 1 ) (erg cm^s" 1 ) 



May 2002 


0±o 


35 +3 3 


list™ 




















bb+brems 


8.0 


X 


10" 


■12 


1.4 


X 


10" 


ii 


1.49 


bb+mek 


0±l 


35+\ 


98±i° 


8.0 


X 


10" 


■12 


1.4 


X 


10" 


11 


1.46 


reflect (bb+mek) 


t)+2 
U -0 


351 1 


3o±;g 


8.3 


X 


10" 


■12 


7.2 


X 


10" 


12 


1.44 


absori(bb+mek) 


o+7 
J -3 


32±* 




1.4 


X 


ltr 


■11 


6.7 


X 


10" 


12 


1.37 


reflect * absori(bb+mek) 


0+ 7 
u -o 


34i 3 


i5i 6 3 


1.5 


X 


10" 


■11 


5.6 


X 


10" 


12 


1.36 


May 2003 


























bb+brems 


1 s+ 4 


9Q+ 4 


cn+40 
O 9 -30 


9 R 


X 


1U 


■11 


1 1 
1.1 


X 


1U 


11 


1.4D 


bb+mek 


2.0 +4 2 


29 + 3 


04 -30 


2.7 


X 


10" 


■11 


1.0 


X 


10" 


11 


1.43 


reflect (bb+mek) 


2.7+1 


29 + 3 


22+1° 


3.4 


X 


10" 


■11 


5.9 


X 


10" 


12 


1.38 


absori(bb+mek) 


5.9± 3 3 


25t 3 


15t 5 3 


8.7 


X 


10" 


■11 


5.7 


X 


10" 


12 


1.39 


reflect * absori(bb+mek) 


5.8l 6 3 


25l 3 3 


13l 3 2 


8.6 


X 


10" 


■11 


5.1 


X 


10" 


12 


1.36 


Oct 2003 


























mek 


n+10 
z -2 




3.5± 3 










4.8 


X 


10" 


14 


0.5 


May 2005 


n+10 
U -0 
























mek 




1.96±i 










4.3 


X 


lQ- 


14 


1.25 



1 The JVh refers to cold absorption applied to all model components. 



temperature of the reprocessed X-ray emission. One par- 
ticul ar example is the ROS AT PSPC observation of AM 
Her (jBeuermann et aLll2008l ), which shows a correlated in- 
crease of the blackbody temperature with increasing soft 
X-ray flux. This result is an indication that the tempera- 
ture is set in response to a locally enhanced accretion rate, 
at least for this given system at certain epoch. The anti- 
correlation of temperature and soft X-ray flux in HU Aqr 
would contradict this picture and could be explained by a 
much smaller soft X-ray emitting area in the intermediate 
state. 

During the intermediate state the main accretion spot 
shifted toward the line connecting both stars to an azimuth 
of ip = 31°. It thereby follows the accretion-rate dependent 
migr ation previously seen in some ROSAT/EUVE point- 
ings (Schwopc ct al. 2001). We did not observe the irregular 
behaviour seen in a few ROSAT/EUVE intermediate state 
observations, where the accretion spot flipped back to the 
high state position. The accretion stream was detected via 
the absorption dip. Compared to the high state its column 
density was reduced by a factor of 10, while its position 
moved by ~ 25° with respect to its high state position, also 
in accord with previous observations. 

The origin of X-rays during the low state of HU Aqr re- 
mains ambiguous. While the luminosity is consistent with 
a coronal emitter, the observed temperatures are better ex- 
plained by an accretion plasma. For the latter interpreta- 
tion we would expect a clear orbital modulation which is 
not observed. 

The (O — C) diagram of the accretion spot eclipse tim- 
ings in HU Aqr reveals large deviations from a constant pe- 
riod on relatively short timescales. It should be noted that 
the observed spot timings, beside a true change of the or- 
bital period, are influenced by the location of the accretion 
spot on the surface of the white dwarf, which is established 
in HU Aqr to depend on the mass accretion rate. Taking 
the observed accretion geometry changes from high to in- 
termediate states, the movement of the accretion spot from 



i> = 45° to tp = 35° (|Schwope et al.ll2001h . we estimate the 
possible, additional shift of the (O — C)-values to be only 2 
sec, being thus small compared to the observed deviations. 

The observed variations in the (O — C)-diagram must 
either be attributed to a secular period variation due to an- 
gular momentum loss or cyclic variations or a combination 
thereof. The situation remains so far ambiguous with both 
a quadratic and sinusoidal ephemeris being equally possi- 
ble and providing a similarly bad or good description of 
the data. A combined quadratic and sinusoidal ephemeris 
results in a significantly better fit. If true, the observed pe- 
riod decrease (P or b = — 7.. — 11 x 10~ 12 )s s~l would be large, 
implying an angular momentum loss of J = 3.7 x 10 35 erg 
neglecting the effects of mass loss. This is in the ran ge of 
standard magnetic braking (jRappaport et all [198 3) . but 
a factor of 30 larger than expected by gravitational ra- 
diation alone. Such a strong angular momentum loss for 
a system below the period gap would clearly violate the 
standard evolutionary scenario of disrupted magnetic brak- 
ing. Notably, a period decrease at a similar level was re- 
cently detected in thr ee other s hort period binaries: DP Leo 
(ISchwope et a"l]|2002l). NN Ser (|Brinkworth et al.ll2006fl and 
OY Car (|Greenhill et al.|[2006T ). Future observations have to 
show if these variations are indeed secular or if their peri- 
odic nature is masked by the yet insufficient baseline. 

iBorges et al.l (|2008l) have compiled a selection of 14 
eclipsing CVs for which a well sampled O — C-diagram has 
been established. All these systems show variations of the 
period, which in the long-run mostly appear to be cyclic 
masking any remaining secular trend. The (O — C)-diagram 
of HU Aqr also reveals clear indication for a cyclic change, 
irrespective of whether a strong secular period evolution 
is assumed or not. The measured modulation periods (13 
or 6.9 yrs) and the amplitudes imply a fractional period 
change AP/P mo j of 3 — 7 x 1 0~ 7 , which is comp arable to 
the values of the four CVs of IBorges et al.1 (|2008[ ) that lie 
below the period gap. The preferred mechanism invoked 
to explain these cyclic variations in CVs and other long- 



8 



R. Schwarz et al.: XMM monitoring of the eclipsing polar HU Aqr 



period close binaries (Algol, RS CVn and W UMa stars) 
are period modulations due to shape cha nges induced by a 
variable magnetic field (|Applegatelll992t ). 

A recent recalculation of App legate's mechanism for 
the short-period pre-CVs NN Ser (|Brinkworth et al J 12006) 
shows that the energy budget of the low-mass secondary is 
one order of magnitude lower than the energy needed to 
drive the observed period variations. The new calculations 
take into account a variable shell mass and the contribu- 
tion of the qu adrupole moment of the core. Following the 
calculations of Brinkw orth et all (|2006[ ) we tested if the en- 
ergy budget of the secondary star in HU Aqr allows for the 
Applegate mechanism at all. Using M2 = 0.2 M , R2 = 0.22 
Rq and Mwd = 0.88 M Q (Vogel et al., in prep.) the change 
in period implied by the sinusoidal terms of the ephemerides 
result in a required energy to drive the Applegate mecha- 
nism which exceeds the total radiant energy of the sec- 
ondary in the case of the quadratic + sinusoidal ephemeris 
by a factor of 3.5 and still requires 80 per cent at min- 
imum in the case of the sinusoidal ephemeris. Neglecting 
any angular momentum loss the quadratic ephemeris also 
requires more than half (~ 70%) of the total energy output 
of the star. We therefore claim that the Applegate mecha- 
nism cannot account for the possible sinusoidal modulation 
of the orbital period. 

As an alternative origin the cyclic period changes could 
also be caused by an unseen third body in which case the 
variation should be strictly periodic. Taking the elements 
of the quadratic plus sinusoidal ephemeris we compute the 
minimum mass of such a body assu ming that its inclination 
is 90° . The observed mass function (Bo rkovits &: Hegeduesl 
1996) is 1.1 x 10 _7 Mq implying a minimum mass of the 
third body of 0.0047M Q (5 Jupiter masses) for an assumed 
total mass of the system of 1.0 M Q . On the otherhand, a 
third body more massive than an object at the hydrogen- 
burning limit would require an inclination < 3° and is 
therefore not very likely. 
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Table 2. Eclipse egress times of HU Aqr. BJED is the barycentrically corrected ephemeris time of the eclipse egress 
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Cycle 


BJED 
_i_9/t /innnn 

i-Z^UUUU 


ABJED 


Instrument 


1319 


9217.4361120 


0.0000115 


MCCP+CA2.2 


1320 


9217.5229220 


0.0000115 


MCCP+CA2.2 


1321 


9217.6097490 


0.0000115 


MCCP+CA2.2 


1322 


9217.6966010 


0.0000231 


ESOlm 


1333 


9218.6516100 


0.0000231 


ESOlm 


1334 


9218.7384390 


0.0000231 


ESOlm 


1367 


9221.603501 


0.0000231 


ESOlm 


1368 


9221.690319 


0.0000231 


ESOlm 


1369 


9221.777148 


0.0000231 


ESOlm 


14087 


10325.959254 


0.0000277 


HST 


14088 


10326.046074 


0.0000277 


HST 


14236 


10338.895480 


0.0000277 


HST 


13620 


10285.414104 


0.0001157 


AIP 


13621 


10285.500878 


0.0001157 


AIP 


13632 


10286.455890 


0.0001157 


AIP 


14115 


10328.390324 


0.0001157 


AIP 


14116 


10328.477162 


0.0001157 


AIP 


14138 


10330.387041 


0.0001157 


AIP 


14139 


10330.473866 


0.0001157 


AIP 


22788 


11081.383711 


0.0001157 


AIP 


27394 


11481.278636 


0.0000231 


AIP 


29955 


11703.625705 


0.0000810 


AIP 


29966 


11704.580704 


0.0000810 


AIP 


35043 


12145.367925 


0.0000463 


AIP 


29955 


11703.625705 


0.0000925 


CA123 


29966 


11704.580704 


0.0000925 


CA123 


31312 


11821.441021 


0.0000115 


STJ 1 


31313 


11821.527841 


0.0000115 


STJ 1 


25892 


11350.874322 


0.0000115 


Optima 


25926 


11353.826219 


0.0000115 


Optima 


25938 


11354.868078 


0.0000115 


Optima 


30276 


11731.495069 


0.0000115 


Optima 


30277 


11731.581892 


0.0000115 


Optima 


35376 


12174.279097 


0.0000115 


Optima 


35377 


12174.365920 


0.0000115 


Optima 


38109 


12411.559191 


0.0000115 


Optima 


39731 


12552.381844 


0.0000115 


Optima 


39742 


12553.336841 


0.0000115 


Optima 


42441 


12787.665041 


0.0000115 


Optima 


42463 


12789.575075 


0.0000115 


Optima 


42464 


12789.661918 


0.0000115 


Optima 


47253 


13205.444710 


0.0000173 


Optima 


47254 


13205.531527 


0.0000173 


Optima 


47300 


13209.525279 


0.0000173 


Optima 


47335 


13212.563996 


0.0000173 


Optima 


48265 


13293.306950 


0.0000173 


Optima 


48288 


13295.303800 


0.0000173 


Optima 


48299 


13296.258845 


0.0000173 


Optima 


48334 


13299.297557 


0.0000173 


Optima 


51032 


13533.539015 


0.0000173 


Optima 


55466 


13918.500708 


0.0000173 


Optima 


55546 


13925.446340 


0.0000173 


Optima 


60096 


14320.479250 


0.0000173 


Optima 


60097 


14320.566050 


0.0000173 


Optima 


38105 


12411.211871 


0.0000578 


XMM OM 


42352 


12779.938034 


0.0000578 


XMM OM 


38107 


12411.385572 


0.0000231 


XMM MOS1 


38107 


12411.385578 


0.0000231 


XMM MOS2 


38108 


12411.472392 


0.0000231 


XMM PN 
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Table 2. continued. 



Cycle 


BJED ABJED 


Instrument 




+2440000 





38133 12413. 
38145 
42395 
50702 
50713 
50714 
50724 
50725 
50737 
59524 
59525 
59558 
59559 



642851 
12414. 
12783. 
13504. 
13505. 
13505. 
13506. 
13506. 
13507. 
14270. 
14270. 
14273. 
14273. 



0.0000115 
684740 0. 
671299 
888294 
843317 
930139 
798342 
885162 
927008 
817989 
904812 
769884 
856704 



ULTRACAM 

0000115 ULTRACAM 

0000115 ULTRACAM 

0000056 ULTRACAM 

0000056 ULTRACAM 

0000056 ULTRACAM 

0000056 ULTRACAM 

0000056 ULTRACAM 

0000056 ULTRACAM 

0000157 ULTRACAM 

0000157 ULTRACAM 

0000157 ULTRACAM 

0000157 ULTRACAM 



(1) From Bridge et all £2002) 



